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The periodic and step-like solutions of the double-Sine-Gordon equation are investigated, 



{^JQi with different initial conditions and for various values of the potential parameter e. We plot 

energy and force diagrams, as functions of the inter-soliton distance for such solutions. This 



allows us to consider our system as an interacting many-body system in 1 + 1 dimension. 
We therefore plot state diagrams (pressure vs. average density) for step-like as well as 
periodic solutions. Step- like solutions are shown to behave similarly to their counterparts 
in the Sine-Gordon system. However, periodic solutions show a fundamentally different 



behavior as the parameter e is increased. We show that two distinct phases of periodic 
solutions exist which exhibit manifestly different behavior. Response functions for these 
phases are shown to behave differently, joining at an apparent phase transition point. 
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^ The Sine-Gordon (SG) equation is a non-linear partial differential equation which appears natu- 

^ I rally in different physical systems in atomic physics electromagnetismy], superconductivity Q], 



systems : 

field theory [4], biophysics [a, I^, Q], and statistical mechanics [ 

The double-Sine-Gordon (DSG) equation which is a generalization of the ordinary SG equation 
has been the focus of much recent investigations. It has been shown to model a variety of systems in 



condensed matter, quantum optics, and particle physics 



the spin dynamics of superfluid ^He 



Id, 



icslsl. 



Condensed-matter applications include 



111, magnetic chains il2i1. commensurate- incommensurate 
Phase ...i«o„s|li ...face *„c.u», domain waUsQ Q aM 

namics in Josephson j unction [l 



In quantum field theory and quantum optics DSG applications include quark confinement 
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FIG. 1: DSG Potential. The dashed curve is for e = 10, the dash-dotted curve is for e = 1 and the sohd 
curve is for e — 0(SG). 

and self-induced transparency [l^ . The internal dynamics of multiple and single DSG soliton 
configurations using molecular dynamics have been studied in[9|. There have also been studies 
about kink anti-kink collision processes for DSG equation [20^. One can also point to the statistical 
mechanical applications 2l|], and perturbation theory for this equation [22]. It should be mentioned 
that the potentials adopted in various studies are not exactly the same. The DSG potential which 
contains a constant and a harmonic term in addition to the self-interaction potential of the ordinary 



SG equation is considered here 
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V{(p) = l-he-cos0-ecos(2(/>). 



(1) 



where e is a constant. This potential is sketched in Fig{T] for e = 0, 1 and 10. Some dynamical 
properties of multiple and single soliton solutions of the DSG system were studied by Burdick et. 
al[9|, using molecular dynamics. 

Following the common terminology in field theory, this potential has absolute degenerate minima 
at <j) = 2mT as the true vacuua, and the metastable, local minima at <j) = (2n + l)7r as the false 



vacuua 



, 23(1 . False vacuua develop if e > 0.25[23]. The harmonic term in this potential can result 
from the Fourier expansion of an arbitrary, periodic potential V{(f>) = Vifp + 2n7r). One does not 



23|. The potential reduces to 



expect the system to remain integrable by adding these extra terms 
the ordinary SG potential in the limit e — > 0. 

In this article, we consider chains of the kink and anti-kink solutions of DSG equation in 1 -|- 1 



3 



dimension. Our goal is to consider a one dimensional chain of many solitons and treat them 
as a many-body interacting system. As will be shown shortly, these solutions do have mutual 
interactions which is a function of their mutual distance. We thus propose to study the DSG 
system as a many-body interacting system at zero temperature in one spatial dimension. 

The structure of our presentation is as follows: in Section [TI] we review some basic properties of 
the DSG system. In Section Hill we study the periodic and step-like solutions of the DSG system, 
highlighting their interaction properties. In Section IIVI we study the macroscopic properties of 
such systems viewed as a many-body interacting system. We close in Section |V] by summarizing 
our results and pointing some directions for future work. 

II. BASIC PROPERTIES OF THE DOUBLE SINE-GORDON SYSTEM 

From a relativistic point of view, the double-Sine-Gordon Lagrangian density can be put in the 
following form: 

Cdsg = ^d''(l)d^<J) - [1 + e - COS./. - ecos(2(^)]. (2) 

From this Lagrangian density, we obtain the following equation: 

= -sin(/)- 2esin(2,/>); (3) 

for the real scalar field (j){x,t) in (1 + 1) dimensions. It can be shown that this equation has the 
following exact, static, single kink (anti-kink) solution[23^: 



(/>(x) = 2 arccos 



sinh v'4e + Ix 



(4) 



V^e + coshA/ie+Tx 

These are plotted in Figl2] for three values of the parameter e(= 0,1,10). Note that the (+) 
sign is for the anti-kink (k), while the (— ) sign is for the kink (k). The solution ^ tends to the 
conventional kink (anti-kink) solution of the SG equation in the limit e — > 0. As the value of e 
grows, the kink develops two parts (or sub- kinks) I23I , see FigjSl 

Exact solutions of DSG system with arbitrary constant coefficients were obtained recently by 
Wang and Lie[24] by F-expansion method which is a generalization of the Jacobi elliptic function 
expansion. It is obvious that solitary and periodic solutions of the SG system are special cases 
of these solutions (with particular choices of the constants of integration). The thermodynamic 
properties of the DSG chain were studied in Condat et. al[25,], in which the polarization precursor 
was determined. 
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FIG. 2: Kink and anti-kink solutions of DSG. The dotted curves are for e = 10, the dash-dotted curves are 
for e = 1 and the dashed curves are for e — 0(SG). 

Using the Noether's theorem and the invariance of the action under the the space-time trans- 
lation -^^ -|-a'^, one easily obtains [3]: 



(5) 



in which g^" = diag{l, —1) is the metric of the (1 -|- 1) dimensional Minkowski spacetime. The 
energy-momentum tensor of the true vacuum ((p = 2mr) vanishes, while for the false vacuum 



(2n -|- l)vr), it is non- vanishing 



(6) 



231, 



As expected, this tensor is that of a perfect fluid, with the equation of state p = —p = 
similar to that of the cosmological constant [27I . Using the virial theorem, the energy density of 
the static kink is given by (see FigJS]): 



nix) = 2V{(\)) = 2[l-h€-cos((/>(rE))-ecos(2</.(x))]. 
As for the SG system, the topological current for the DSG system is given by[28^: 

where e^^ is the totally antisymmetric tensor in two dimensions. The current density 
served: 

1 



(7) 



(8) 



IS con- 



Mil 



2tt 



d^d^d) = 0. 



(9) 
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FIG. 3: DSG energy density for the kink solution. The dashed curve is for e ~ 10, the dash-dotted curve is 
for e = 1 and the solid curve is for e — 0(SG). It can be seen that sub- kinks develop for large enough values 
of e (e > 0.25) 

since e^'^ is antisymmetric while df^du is symmetric in the indices. The total topological charge of 
a localized solution is given by: 

jOrfx = — / €°^di(l)dx 

-oo '^'^ J—oo 

Considering the boundary conditions of a localized solution, the DSG equation leads to: 

(/)(-|-oo) = 2mr, (j){—oo) = 2rmr; Q = n — m. (11) 
Localized, static solutions have n — m = ib27r, and therefore 

Q = ±l for kink (anti-kink). (12) 



As the value of e grows, sub-kinks form, which have half- integral topological charges|23l|. 

The stability of kink solutions of the DSG system is guaranteed by topological reasons. The 
kink is the lowest energy configuration with the prescribed boundary conditions and the topological 
charge given by Eq. (fT2]l . 



III. PERIODIC AND STEP-LIKE SOLUTIONS 



In this section, we classify and study two types of static solutions of the DSG equation, for 
which we have: 
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I 1 

In order to proceed, we use a Runge-Kutta method [29] to integrate Eq. (I13p numerically. We 
choose (j) = TT and various values of ^ as our initial conditions at x = 0. We study this system for 
various potential parameters e. We choose e = 1 as a typical value but also compare results with 
e = 0.1 (SG-like behavior) as well as e = 10 for decidedly different behavior. It is easy to see that 
solutions are characterized by an (x-independent) constant P: 

^-Ufi-vi^y. (14) 



dx 

Note that this is different from energy density 



(15) 



1 fdsV 



which changes with position x. We will interpret P as "pressure" (or tension in this Id case) for 
reasons to be explained in the next section. 




FIG. 4: Step-like chain of DSG solitons for e = 1 and P = 0.42. The inset shows the cyclic nature of these 
solutions in more details. 
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FIG. 5: Periodic chain of DSG solitons for e = 1 and P ~ —1.8848. The inset shows the cychc nature of 
these solutions in more details. 

Two types of solutions emerge under such conditions: Step-like solutions and periodic solutions 
(see FigH] and [5|). Step-like solutions are characterized by P > and periodic solutions are 
characterized for —2 < P < 0. Note that this classification is independent of e and therefore holds 
true also for SG equation[30]. As can be seen from the relevant figures the step-like solutions are a 
sequence of kinks (kkk...) or anti-kinks (kkk...), and the periodic solutions are a sequence of kink, 
anti-kink {kkkk...) solutions. A simple phase diagram of Eq. (13) is shown in Fig. 6. 

Using Eq. ()16p . one can characterize these solutions via an energy function. We take a given 
static solution and compute the total energy per "soliton" i.e. integral of Eq. (fT6l) over one period. 
For each static solution we also calculate the "inter-soliton" distance, L, as a distance between two 
successive (similar) peaks of energy density. We can therefore plot energy per soliton as a function 
of inter-soliton distance for various values of e and P. Since it is our goal to consider these systems 
as a many-body system, we use energy as a function of distance E(L) and compute the "inter- 
soliton force" F = Fi gure [7] shows such results for the step-like solutions. No significant 
change of behavior is seen as e is varied for these type of solutions. We note that in both cases 
(SG and DSG) at short distances the energy rises abruptly reminiscent of hard core potentials. 
It quickly falls as distance L is increased reaching an asymptotic value for large L. Accordingly, 
the "inter-soliton force" F (insets) is repulsive and falls to zero quickly as distance is increased. 



-6i- 
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FIG. 6: A typical phase diagram for solutions of Eq. (13) with e = 1.0. The two different solutions are 
clearly marked on the diagram. 



indicating a short-range interaction. 

We next consider the periodic solutions. Here, as we will show, e changes the behavior of these 
solutions in a significant way. Figure [8] shows the results for the periodic solutions for e = 
(SG). The behavior is the reverse of the step-like solutions. The energy decreases sharply at short 
distances and increases quickly to a constant value at large L. Therefore, the inter-soliton force is 
attractive and falls off to zero at large distances (inset). What happens as e is changed in this case? 
The results for the energy diagram is shown in FigJOj As is seen from the figure, the introduction 
of e leads to the emergence and growth of a new branch of solutions. It is important to distinguish 
between these two types of periodic solutions. The inter-soliton distance L changes with P which 
is a function of initial conditions. In the lower branch (the SG-like branch), as L becomes large, 
the energy per soliton tends to its single-soliton value, thus leading to an independent soliton 
system. However, in the upper branch, as L become large the sub-kinks also separate out with 
a false vacuum energy developing between them. Here, as L grows large, the total energy in the 
false vaccua increases linearly in the space between the sub-kinks, which is indicated by the linear 
dependence of energy in Figjolsi]. To help visualize these different solutions, we show in FigllOl a 
periodic solution from each branch, having the same L but different E. The inter-soliton force for 
the periodic solution is still attractive since in either branch the energy is an increasing function of 
distance. However, in the lower (SG like) branch the force falls off to zero at large distances leading 
to a non-interacting soliton picture. Note that this is the case in all the solutions considered thus 



9 




FIG. 7: Energy per soliton diagram for step-like chain of (a) SG solitons and (b) DSG solitons for e — 10. 
The behavior of the system dose not change significantly as e is changes. The inset shows the derivative 
which we interpret as the (repulsive) force between the solitons. 

far. But, in the new upper branch, the attractive force increases its (absolute) values to —2 as L 
is increased (Fig. 10b). This is due to the above-mentioned linear behavior of energy in the upper 
branch, caused by separation of sub-kinks with increasing L. A typical force diagram for periodic 
solution of DSG system with e = 1 is shown in Figllll 



IV. PERIODIC AND STEP-LIKE SOLUTIONS AS A MANY-BODY SYSTEM 



In order to gain a better understanding of the soliton interpretation and mutual interactions 
between solitons, one may also study state diagrams which depict pressure P versus average density 
p, where P and p are given by 32] : 

2 



(17) 



It should be noted that the energy-momentum tensor of the system (like many other relativistic 
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FIG. 8: Energy per soliton diagram for periodic chain of the SG sohtons. The inset shows the derivative 
which we interpret as the (attractive) force between the sohtons. 
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FIG. 9: Energy per sohton diagram for periodic chain of DSG sohtons (a) for e = 0.1, (b) for e = 1 and 
(c) for e = 10. Note the emergence of a new (upper) branch as e is increased. The lower branch reaches a 
constant value for large L, but the upper branch increases linearly with slope 2. 

continuous media) takes the form 



TO 
Tl 



(19) 



reminiscent of the energy-momentum tensor of a perfect fluid in 1 + 1 dimensions given by 
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FIG. 10: The slope diagram for the periodic chain of DSG sohtons for e = 1 , with P = —2.14 x 10^^ for 
the curve with larger amplitude (solid curve) and with P = —1.9996 for the curve with smaller amplitude 
(dash-dotted curve). The two solutions have the same inter-soliton distance (L=7.798) but different energy 
per soliton. 




FIG. 11: Attractive force for: (a) the lower branch and (b) the upper branch of periodic chain of DSG 
solitons for e = 1. 



(Tn = I ° I (20) 

-P 

Note that pressure is the same as our first integral P (Eq. (|14|) ) and density p is in fact the same 
as energy (Hamiltonian) density (Eq. (jl6p ). Here we use the generic word "pressure" as a thermo- 
dynamic conjugate variable to the size of the system. If the system is three dimensional, pressure 
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has its usual meaning. However, for one dimensional systems like the ones we are considering here 
"tension" is a more appropriate name for our variable, {rdl instead of —Pdv as mechanical energy 
or work). With this in mind, we use the word pressure or tension interchangably. 

In Fig ll2l we show a state diagram for SG as well as DSG for e = 10 for the step-like solutions. 
We note that, as before, the behavior is essentially the same. The pressure increases smoothly 
from zero and tends to a linear regime with increasing density. This behavior is consistent with 
the repulsive force diagram (see Figl7|) which shows that the inter-soliton force becomes stronger 
at short distances which indicates that pressure increases with increasing density. 



P 3 




P 3 




FIG. 12: Equation of state diagram for step-like chain of (a) SG solitons and (b) DSG solitons for e — 10. 




FIG. 13: Equation of state diagram for periodic chain of SG solitons. 
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FIG. 14: Equation of state diagram for periodic chain of DSG solitons for e = 1. For examples of solutions 
from each branch see Fig. 10. 

As the energy and force diagrams indicate, the state diagrams for periodic chains are more com- 
phcated (and therefore more interesting) than the step-like solutions. In Fig llSl we show the state 
diagram for a periodic chain of SG solitons. Here, with increasing density, the pressure decreases 
from zero. Clearly this negative pressure is a result of attractive force (see Figj8]) among periodic 
kink-antikink SG solitons. Note that as density increases and inter-soliton distance decreases, the 
attractive force becomes stronger thus lowering the pressure or increasing the tension in this one 
dimensional chain. The negative value of pressure is simply related to the negative inter-soliton 
force. An interesting feature of this diagram occurs at p ~ 2 which is the maximum density. Here, 
P = —2 and ^ = 0, i.e. the "fluid" becomes incompressible. This simply corresponds to an 
exceedingly large (attractive) force at short distances. We also note that since ^ is negative, 
the term "inexpandable fluid" is perhaps more appropriate than incompressible fluid. Note that 
the end point {P, p) = (—2,2) exactly corresponds to the vacuum equation of state P = —p, 
coming from Eq.Q. Note, however, that this point is highly unstable in the SG system, since it 
corresponds to a local maximum in FiglH 

Next, we study the state diagrams for a periodic DSG system. As shown above, the energy 
diagram splits into two (upper and lower) branches (see FigE]). Figure E] shows a typical state 
diagram for DSG periodic chain. In fact, the state diagram also contains two (connected) branches. 
The upper half which is similar to the previous state diagram for the simple SG system, as well as 
a new, lower half which returns the same values of density but at lower pressure values or larger 
tension. We note that the upper half corresponds to the lower (SG-like) energy branch, whereas 
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the lower half corresponds to the new upper energy branch of DSG system. Again, it is easy to 
understand this state diagram in terms of its corresponding force diagram. Here, in the new lower 
branch, increasing the density (and thus reducing L) causes a weaker attractive force (see Fig JTT]) 
which in turn causes an increase in the pressure. 

The shape of the state diagram, Fig |14l with double-valuedness of pressure is reminiscent of 
Van der Waals fluid. There, the multiple-valuedness of density as a function of pressure with 
a corresponding change in the sign of compressibility X = ^ signals the onset of a liquid-gas 



transition as temperature is lowered 



331]. Our system is obviously very different from that described 



by the Van der Waals equation of state. However, our system exhibits two distinct solutions, one 
with a positive and the other with negative being joined at a maximum density pmax = 3.385. 
This point could be thought of as a phase transition point. The transition point is the point at which 
the "fluid" becomes "inexpandable" , which separates regimes of different signs of compressibility. 
In other words, the transition point is the point at which a stable solution (x > 0) gives way to an 
(thermodynamically) unstable solution {x < 0). This can be better seen in FigllSI where we plot 
^ vs. P in order to show these two distinct phases. 

Energy per soliton (FigJO]), and equation of state diagrams (Fig ll4p for the periodic chain clearly 
show two branches. These two branches correspond to two different phases of the chain. The two 
phases differ in the following respect: The upper branch of Figll4l corresponds to a chain ...kkkkk... 
with solitons separated by a region of true vacuum (zero energy density). This is why at large 
inter-soliton distances the force (and thus tension) vanishes for this branch. The lower branch, on 
the other hand, corresponds to pairs of sub-kinks (with vanishing topological charge) separated by 
regions of false vacuum (non- vanishing energy density). This explains the asymptotic value —2 at 
low densities (large inter-soliton distances, see figjll(b)). It is interesting to note that the energy 
stored in the false vacuum increases linearly with the inter-soliton distance (at large distances) 
which leads to a constant force (and thus tension). This interesting behavior is reminiscent of the 
confinement phenomenon in hadronic physics [ 



V. SUMMARY AND CONCLUSIONS 



We have investigated the periodic and step-like solutions of the double-Sine-Gordon equation in 
this paper. Runge-Kutta algorithm was employed to integrate the static, second order ODE, with 
adjustable initial conditions {(p = n and various ^ at a; = 0). The solutions fall into two categories: 
periodic, and step-like. Step-like solutions are a sequence of kinks {...kkk...) or anti-kinks (...kkk...). 
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FIG. 15: i versus P, showing clearly two phases with positive and negative values of i. 

Periodic solutions are a sequence of kink, anti-kink (...kkkk...). For initial conditions considered 
here, the constant "pressure" P characterizes these solutions with P > for step-like solutions and 
—2 < P < for periodic solutions. We also characterize these solutions using their energy density 
and therefore calculate energy as a function of distance. We are therefore able to consider an 
"interaction energy" as a function of distance. In this regard, we find that the step-like solutions of 
DSG equation are similar to their SG counterparts, i.e. the behavior of the system does not change 
substantially with parameter e. However, the behavior of periodic solutions depends crucially on 
e. In fact, we observe the emergence of an extra branch of solutions as e is increased. 

Using the concept of "interaction energy" we are able to consider our system as an interacting, 
many-body system on a one dimensional chain at zero temperature. Using standard definitions 
of pressure and density, we calculate the equation of state for such a system. In this regard, 
the step-like solutions show standard behavior. However, the periodic solutions show unusual 
behavior. First, due to the nature of attractive force between them, they exhibit negative pressure. 
More interestingly, they exhibit a transition from a region of positive compressibility to a region of 
negative compressibility signifying an intrinsic instability common in many thermodynamic systems 
exhibiting phase transition. This is of particular interest since most one dimensional interacting 
systems do not exhibit thermodynamic phase transition. 

Many interesting questions arise: Is there an experimental realizations for such a system? We 
believe that given the prevalence of DSG equation in many experimentally realizable situations 
outlined in the introduction, there is hope to study such systems experimentally and test our 
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results. 

Furthermore, other generahzation are possible. For example, we have looked at potentials of 
the form V{(/)) = 1 + e — cos(0) — ecos(n0) with n integer larger than 2, i.e. the multiple Sine- 
Gordon equation. Here, instead of two branches we obtain multiple branches. Another avenue of 
investigation is the study of inhomogeneity (e.g. e T^constant) in these system. Finally, the role of 
temperature and fluctuations might be of interest in these systems. Here, the dynamical stability 
of the solutions we have studied here become an important issue. These issues are currently under 
investigation and we intend to report our results in forthcoming publications. 
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